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(NASA-CR-18898 3) MULT I-DI SC I PL I NARY 
OPTIMIZATION OF AEROSERVOELASTIC SYSTEMS 


N92-I1377 


The purpose of the research project "Multi-Disciplinary Optimization of 
Aeroservoelastic Systems" in its second year (from October 1990 to September 
1991) was to continue the development of new methods for efficient aeroservo- 
elastic analysis and optimization. The main target was to develop a method for 
investigating large structural variations using a single set of modal 
coordinates. This task has been accomplished by basing the structural modal 
coordinates on normal modes calculated with a set of fictitious masses loading 
the locations of anticipated structural changes. 

The first part of this work, "Modal Coordinates for Aeroelastic Analysis 
with Large Local Structural Variations", has been presented at the European 
Forum on Aeroelasticity and Structural Dynamics, Aachen, Germany, June 1991. The 
paper was co-authored by Carol D. Wieseman of NASA Langley. The abstract is 
given in Appendix A. The paper has also been submitted for publication in the 
AIAA Journal of Aircraft. The second part of this work, "Time Simulation of 
Flutter with Large Stiffness Changes", has been proposed, together with Carol 
Wieseman, for presentation at the 33rd Structures, Structural Dynamics and 
Materials Conference, April 1992. The submitted extended abstract is given in 
Appendix B. 

Another work on sensitivity derivatives for residualized aeroservoelastic 
optimum-design models has been performed together with Israel Herszberg of RMIT 
Victoria University of Technology, Australia, while he was on sabbatical at the 
Technion. The abstract of the paper on this work that has also been presented in 
the Aachen Conference mentioned above is given in Appendix C. 

Two other subjects that I worked on together with graduate students are 
aeroservoelastic optimization with continuous gust response constraints (with 
Arie Zole) and dynamic response to impulsive excitation (with Eyal Presante). 
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The target in both cases is to develop efficient numerical procedures for 
inclusion in the multi-disciplinary optimization scheme. It is anticipated that 
the first subject will be incorporated in the FASER optimization program by 
June, 1992, and the second subject before the end of 1992. 

A preliminary investigation has been started, in cooperation with Shlomo 
Taason of ICASE on the inclusion of static aeroelastic effect in CFD shape 
optimization. Initial guidelines for 2-D optimization are given in Appendix D. 
The basic principles will be applied to 3-D optimization at a later stage. 
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MODAL COORDINATES FOR AEROELASTIC ANALYSIS WITH LARGE LOCAL 
STRUCTURAL VARIATIONS 

M. Karpel C.D. Wieseman 

Faculty of Aerospace Engineering NASA Langley Research Center 

Technion-Israel Institute of Technology Hampton, Virginia 

The investigation of the effects of local structural property variations on flutter and aeroelastic 
response is an important pha^e of the design of flight vehicles. Common aeroelastic analyses start 
with the calculation of a limited set of normal vibration modes’ which serve as generalized 
coordinates. Repeated calculations of the normal modes and the associated generalized aerodynamic 
force coefficients every time a structural property changes may be impractical. An additional 
problem is how to deal with the changing coordinates during time-domain response simulations. These 
problems can be, resolved by performing the entire analysis with a constant set of modal coordinates. 
Various applications used the vibration modes of a nominal structure as generalized coordinates and 
account for structural changes by introducing stiffness, mass and damping coupling terms. However, 
the application of this approach to cases of large structural variations may either yield inaccurate 
results or require a very large number of modes. 

The approach taken in this work is to calculate the constant set of modal coordinates with the free- 
free structure loaded with fictitious masses at the vicinity of the varying stiffness. The 
application of relatively large fictitious masses causes significant local deformations in the low 
frequency modes. This facilitates high accuracy usage of a small number of modal coordinates over a 
wide range of local stiffness variations. The fictitious masses are removed in the analysis by 
appropriate mass coupling terms. Damping coupling terms are introduced to yield the desired modal 
damping values at air off conditions. 

The numerical examples deal with a realistic aeroelastic system of an aircraft with a tip store. 
The wing-store pitch stiffness properties are subject to change by a factor of 30, which has a 
drastic effect on the flutter characteristics. It is first shown that without fictitious masses the 
transition-by-coupling from one configuration to the other causes significant errors of critical 
natural frequencies even with 60 modes taken into account. The introduction of a single fictitious 
mass yields very accurate frequencies with as few as 10 modes taken into account. A parametric study 
shows that the results are insensitive to the value of the fictitious mass provided that it is above 
a certain value but not large enough to cause ill-conditioning. The^ state-space aeroelastic 
equations of motion are formulated using minimum-state rational approximations of the unsteady 
aerodynamic forces. The differences between flutter velocities and frequencies obtained by coupling, 
and those obtained by direct solutions are less than 3%. 


TIME SIMULATION OF FLUTTER WITH LARGE STIFFNESS CHANGES 

Mordechay Karpel 

Technion-Israel Institute of Technology 
Haifa, Israel 

Carol D. Wieseman 
NASA Langley Research Center 
Hampton, VA 23665-5225 

INTRODUCTION 

The common approach for formulating the equations of motion of aeroelastic systems 
starts with normal modes analysis of the structural model. Unsteady aerodynamic force 
coefficient matrices are then calculated at various reduced frequency values to account for 
the aerodynamic coupling between these modes while undergoing oscillatory motion. Classic 
frequency-domain aeroelastic analysis methods use the modai structural properties and the 
tabulated aerodynamic matrices for calculating the flutter conditions at which the aeroelastic 
system matrix is neutrally stable. The basic assumption of the modal approach is that the 
structural deflections of the aeroelastic system are linear combinations of a limited set of low 
frequency vibration modes. , ... , . c . 

Time-domain aeroelastic modeling techniques, which cast the equations of motion in a 
state-space, time-invariant form, call for the approximation of the aerodynamic matrices by 
rational functions in the Laplace domain. The order of the resulting state-space model is a 
function of the number of selected modes, the number of aerodynamic approximation roots, 
and the approximation formula. The main considerations in constructing the model are its 
size (which affects the efficiency of the subsequent analyses), its accuracy, and the model 
construction efforts. 

Various computational schemes such as structural optimization, parametric studies, the 
investigation of damage effects, and structural changes during dynamic response, require 


repeated construction of the model for numerous structural variations. Repeated calcula- 
tion of the normal modes and the associated aerodynamic matrices every time the structure 
changes is often impractical in these cases. A more practical approach is to introduce struc- 
tural changes without changing the modal coordinates. The validity of this approach depends 
on the structural information contained in the modal coordinates and the magnitude of the 
structural changes. The number of required modes increases with the magnitude of the 
allowable structural move limits. Keeping the modal coordinates unchanged is even more 
important in time simulation of aeroelastic response during which structural changes occur. 
The occurrences of structural changes define the time segments between which the model 
changes. The end conditions of one segment are the initial conditions of the following one. 
One can argue that these transition-point conditions can be transformed to a set of new co- 
ordinates. But this transformation is adequate only if the new coordinates can be expressed 
as a linear combination of the old ones. However, if this is the case, there is also no need to 
change the modal coordinates in the first place. 

Various structural optimization procedures, like that of Ref. 1, demonstrated that mod- 
erate structural changes can be accommodated .without, changing the modal coordinates. In 
order to accommodate large structural changes, the modes that serve as generalized coordi- 
nates must contain significant distortions in the vicinity of the changes. Taking into account 
more modes supports this purpose but may result in an excessively large aeroelastic model. 
A method for accommodating large structural changes at a small number of structural lo- 
cations, without significantly increasing the model size has been presented in Ref. 2. The 
procedure starts with calculating a set of low frequency vibration modes with the structure 
loaded with large fictitious masses at the location of anticipated structural changes. The 
fictitious masses cause the vibration modes to contain the local deformations required for an 
adequate accommodation of large structural changes. 



The purpose of the proposed paper is to outline the process of using fictitious-mass modes 
to generate efficient fixed-coordinate time-domain aeroelastic models for dynamic response 
and loads analyses, and to demonstrate the usage of these models for time simulation of 
flutter during which large local stifFness changes occur. 

SELECTED RESULTS 

The numerical example consists of a mathematical model of the AFW wind-tunnel model 
tested at the NASA Langley Research Center. An external store is connected to the tip of 
the wing through a decoupling mechanism designed to decouple the pitch motion of the store 
from that of the wing when flutter occurs. With ’’stiff” or ’’coupled” pitch connection, the 
Mach 0.9 antisymmetric flutter dynamic pressure is about 1.9 psi and the flutter frequency 
is about 12 Hz. When the decoupling mechanism is activated, the pitch connection stifFness 
is reduced by 96.5% (’’soft” or ’’decoupled”) and the flutter conditions change to 2.9 psi, 31 
Hz, namely a drastic change in both flutter dynamic pressure and the flutter mechanism. 

The fictitious-mass finite-element (NASTRAN) model is with the soft pitch connection 
and with a fictitious pitch inertia of 3 lb-in-sec 2 , twice that of the tip store, loading the wing 
end of the pitch spring. A set of 14 low frequency vibration fictitious-mass (FM) modes 
is used to generate the basic stiff model. A comparison between the first seven natural 
frequencies obtained directly from NASTRAN and those calculated from the 14 FM modes 
and from 14 soft modes without fictitious masses is given in Table 1. It is clear that the FM 
model produces accurate frequencies while the flutter critical store pitch frequency obtaibed 
from the soft model without fictitious mass is 18% higher that the correct one. 

The 14 FM modes were used to generate a time domain aeroelastic. model for the stiff 
case using Minimum-State rational approximation of the unsteady aerodynamics 3 with 8 
aerodynamic augmenting states. The transition to the soft condition is performed by simply 
introducing the appropriate stifFness coupling terms and damping coupling terms (to yield an 



effective diagonal modal damping matrix). The modal coordinates, as well as the associated 
aerodynamic terms, remain unchanged. A comparison between flutter conditions calculated 
by the frequency domain p-k method and those obtained from the root-locus analysis of the 
state-space models is given in Table 2. While the direct sitiff and soft cases required full 
analyses starting from separate finite-element models, the FM results were obtained from 
models that are based on the same modal coordinates and aerodynamic coefficients. 

The adequacy of the FM model for time simulation is checked by comparing responses to 
initial conditions with those obtained by the direct models. The dynamic pressure at which 
responses are calculated is 2.5% higher than the stiff flutter dynamic pressure. The vertical 
acceleration responses of the foreward and rear tip-store points are given in figure 1 for the 
stiff case and in figure 2 for the soft case. The development of flutter at this dynamic pressure 
and the decay due to activation of the decoupling mechanism are simulated in figure 3. The 
sransition from stiff to soft is performed by introducing the appropriate coupling terms 
between the the structural states of the FM model. The only output equations that need to 
be changed in the transition are those involve acceleration response. It can be observed that 
the system stabilizes very fast but there may be significant overshoots in the transition. 

CONCLUDING REMARKS 

Normal modes calculated with fictitious masses at selected structural locations form a 
set of generalized coordinates with which aeroelastic systems can be analyzed over a wide 
range of stiffness changes in the vicinity of these location. The values of the fictitious 
masses can be arbitrarily chosen from a wide range of masses that cause significant local 
deformations in the low frequency modes, without causing numerical ill conditioning. The 
physically weighted minimum-state rational approximation of the unsteady aerodynamic 
force coefficient matrices is performed once for all the stiffness variations. The size of the 
resulting time-domain model is similar to those required for the anatyses of each stiffness 



case separately. The model facilitates an efficient simulation of aeroelastic time response 
during which large stiffness changes occur. A numerical application demonstrated a flutter 
divergence followed by rapid convergence due to the activation of a decoupling mechanism 
which reduces the stiffness of a critical element by 96.5%. The example demonstrates that 

^ , j • , 

the activation of the decoupler may cause a transient overshoot before the system stabilizes. 
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FIGURES 

1. Time response to initial condition, stiff store connection. 

2. Time response to initial condition, soft store connection. 

3. Time simulation of flutter with stiff-to-soft transition. 


Table 1. Natural frequencies of the stiff model 


Number 

direct 

from NASTRAN 

from soft 
with FM 

from soft 
without FM 

Description 

i 

0 

0 

0 

rigid-body roll 

2 

7.023 

7.024 

7.032 

1st fuselage bending 

3 

7.856 . 

7.864 

7.980 

1st wing bending 

4 

13.069 

13.179 

15.41 

store pitch 

5 

16.161 

16.134 

16.43 

2nd fuselage bending 

6 

27.408 

27.410 

27.47 

3rd wing bending 

7 

38.271 

38.300 

38.76 

1st wing torsion 
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Figure 2. Time response to initial condition, soft store connection 
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Text ► Analytical expressions for the sensitivity derivatives of flutter dynamic pressure and control 

1 ^-spaced design criteria, with respect to structural and control design variables are extended to deal with 
residualized aeroservoelastic mcdels. 

The model considered is a first-order time-domain feedback system subject to model size reduction by 
means of residualization. The structure is represented by a number of structural vibration modes, 
the aerodynamics by rational approximation functions and the control system couples control surface 
motion to the response of motion sensors located on the structure. 

The dynamic residualization technique and the minimum-state approximation method for aerodynamic 
modelling have been demonstrated to yield a relative low order, high accuracy model which may be 
used for on-line design studies where fast response time is important The extended formulation and 
new techniques of this work make the calculations more efficient and further enhance on-line design 
studies. 

To date the optimal design process makes use of the reduced order system ter calculate flutter 
dynamic pressures and gain margins, but reverts to the full model for the calculation of stability 
derivatives. The calculation of sensitivity derivatives depends on the special topography of the 
system matrices which is not pieserved by the residualization process. A large proportion of the 
computation time is spent in calculating these derivatives, particularly if there is a large number 
of design variables. The techniques presented in 1 this work allow the calculation of sensitivity 
derivatives for the reduced order system. The consequent reduction in computation time for the 

stability derivatives is of the order of the square of the ratio of the numbers of states. 

i - ' • • v ; 

A realistic numerical example of the optimal design of an aeroservoelastic system comprising a 

Jcarbon-composite wing with four control surfaces is used to demonstrate the method and to compare 
~ ‘ . jthe. results with those obtained by using sensitivity derivatives computed from the full model. The 

the-'' weight- . with' respect to structtirai" ; arid’ 'control ? variables] 
subject to aeroelastic, control and design, constraints.. * ! 

jThe use of sensitivity derivatives calculated for the reduced system decreases significantly the 

time taken for. the optimization and consequently enhances the optimization method as an on-line 

* ! i 

design tool. j 
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M. Karpel, August 19, 1991 

Aeroelastic Effects in 2-D Shape Optimization 

A given airfoil is subjected to 2-D flow as shown in Figure 1. The "nominal” (or "rigid”) 
angle of attack and trailing-edge control surface deflection are cvo and 8 0 . The x axis is parallel 
to the nominal chord line which connects the leading and trailing edges when a = ao and 
8 = 8 0 = 0. The actual angle of attack (a) is related to the nominal one by 


a = a 0 + a e (1) 

where a e is the elastic deformation of a pitch spring of stiffness k a located at the elastic axis 
at x = cx e , z = cz e and connecting the airfoil to the a setup mechnism (which commands 
oro). The pitch elastic equilibrium equation is 

M = k a a e (2) 

where M is the aerodynamic pitching moment about the elastic axis. The actual control 
surface deflection (5) is related to the nominal one by 

8 — 8o -f- 8 e (3) 

where 8 e is the elastic deformation of the actuator spring k$ representin the actuator stiffness 
and located at x = cz/,, z = czh. The associated equilibrium equation is 

H = k s 8 e (4) 

where H is the aerodynamic hinge moment about the hinge line located at x = ex h . 

The aerodynamic moments are 

M — q f C Pu [x — cx e + z u /(z u — cz e )]dx — q f C pi [x — cx e — z\i(zi — cz e )]dx (5) 
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and 

* 

H = q f C Pu [x — cx h + zj{z u — cz h ))dx - q f C pi \x — cxh — Z\l{z\ ~ cz h ))dx (6) 

Jcx h 

where q = | pV 7 is the dynamic pressure, C p is the surface pressure coefficient, and subscripts 
u and l relate to the upper and lower surfaces respectively. z u and z\ are the upper and lower 
surface coordinates. 

It is assumed that the optimization shape functions include oro, a e , <5 0 and S e as design 
variables, cto is usually utilized by rotating the x-z axes relative to the velocity direction. a e 
defines the shape function 

f o r X j C ( 

The elastic axis and hinge line z locations are defined after geometry changes by the average 
of the local z u and z\ values. The shape function associated with both and <$ e -is 


fs = 


0 


if x < cxh 


(S) 


Xh — xjc otherwise 

The shape functions of Eqs. (7) and (8) are used to define both upper and lower surfaces. 

The elastic parameters a e and 8 e should be treated as dependent variables tuned to satisfy 
the equality constraints of Eqs. (2) and (4). The CFD calculations in each iteration result 
in M and H, and their sensitivity derivatives w.r.t. all design variables including a and S. 
To estimate the changes in a e and 8 e required to satisfy these constraints, we can use 


M 

H 


+ qc 


Cmc C Mt 
Ch 0 Ch 6 
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k a 0 
0. k s 


a e + A a e 

6 e + Ad e 


(9) 


where Cm — M /qc 2 and Ch = H jqc 2 which yields 
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k a 0 

0 ks 


- qc 


Cmc Cm 6 
Cif a Ch s 


({"}-[M]{?}) < io > 


At the end of the iteration, all objective functions and their sensitivity derivatives are 
modified to reflect the effects of the constrain equations, as done for other equality con- 


straints. 



It should be noted that if the optimization is performed for one set of flow conditions 
and required lift, it can be performed without aeroelastic effects. After the optimal ao and 
5 0 are found, M and H are calculated by Eqs. (5) and (6), a e and 8 e are then calculated 
by Eqs. (2) and (4), which yields the actual a and <5 of Eqs. (1) and (3). A multipoint 
optimization, however, and off-optimum calculations, require the inclusion of aeroelastic 
effects. It is suggested to perform an optimization case with a rigid airfoil first (with a and 6 
serving as design variables among others). The elastic deformations can then be calculated 
as described above. Another’ optimization can now be performed for the elastic airfoil when 
the optimal result is known a-priori. This will help studying potential problems associated 
with introducing aeroelastic effects. The formulation given above can be extended to 3D 
optimization by replacing the a e and 8 e parameters by elastic modal deflections that will 


also serve as wing shape functions. 





